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I study a possible effect of momentum deposition from many hard partons traversing the hot and 
dense region produced early in nuclear collisions at the LHC. The expected number of such hard 
partons is large. It is argued that the induced diffusion wakes which carry the momentum de- 
posited to the medium may interact. Due to azimuthally asymmetric geometry in non-central 
collisions, this may lead to small preference of the collective flow in the direction of the impact 
parameter. As a result, a small contribution to the azimuthal asymmetry of hadronic spectra is 
obtained. It may be important to take even such a small contribution into account if quantitative 
conclusions about early thermalisation and low viscosity are to be made, based on the measure- 
ments of the elliptic flow parameter V2. 
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1. Introduction: elliptic flow 

Elliptic flow is the name for an azimuthal asymmetry of hadron production in non-central 
relativistic heavy ion collisions. The name is suggested by the most common interpretation: the 
bulk matter excited to high energy densities expands in the directions transverse to the beam due 
to pressure gradients. In non-central collisions, the overlap region of the two nuclei is not sym- 
metric in azimuthal angle, but rather almond-shaped with the shorter size in the direction of the 
impact parameter. As the pressure gradients are larger where the size of the fireball is shorter, 
stronger transverse expansion will be generated in the direction of the impact parameter. Through 
the Doppler effect, more particles and flatter spectra are then emitted in the direction of stronger 
expansion and this leads to the observed azimuthal asymmetry of the spectra. Hence the name 
elliptic flow. 

Elliptic flow is measured in terms of the second order Fourier coefficient of the azimuthal 
distribution of produced hadrons, V2. This is introduced as 
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where <pR is the azimuthal angle of the reaction plane, defined by the beam axis and the impact 
parameter. Note that other than even cosine terms vanish in symmetric collisions at midrapidity for 
symmetry reasons. 

Observed at RHIC, the elliptic flow is rather large. In hydrodynamic simulations, it turns out 
that such a large flow asymmetry can be achieved only if very fast thermalisation is assumed [Q, ||, 
H|] and the shear viscosity is extremely low [Q]. These are important conclusions from RHIC, though 
their full theoretical understanding is lacking yet. However, in order to ensure the reliability of these 
conclusions one should understand the interpretation of the elliptic flow very well. Particularly, all 
effects which cause elliptic flow in addition to the expansion due to pressure gradients should be 
considered. In this contribution I examine one such possible effect. 



2. Introduction: quenching of hard partons 

In nuclear collisions at highest energy, such as those at the LHC, but to some extent also 
at RHIC, hard scatterings between incident partons happen frequently. They normally lead to 
production of jets or minijets. However, rather few of these jets are indeed observed. The majority 
of the leading hard partons loose energy totally in the surrounding strongly interacting medium. On 
the other hand, this also means that the energy and momentum are transferred to the bulk matter 
and it is appropriate to ask about how the matter responds. The response of the bulk is currently 
widely discussed in the literature []5| |6[ f^, ||, 0]. It is important to note, however, that any such 
response is primarily correlated with the direction of the inducing jet and not with the direction of 
the reaction plane. A possible correlation with the reaction plane may be due to different energy 
loss in various direction of the initial jet. 

At the LHC, we expect many such (mini)jets to be present in a single collision. What is 
the effect on the bulk medium of momentum deposition of all of them? Since the directions of 
their original velocities are distributed isotropically, the simplest expectation would be that after 
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Figure 1: Illustration of the probability that two streams could meet. Left: two streams flowing in the out- 
of-plane direction have a better chance to meet. Right: as the fireball is elongated out of the reaction plane, 
two streams which flow in the in-plane direction have more space to pass each other. 

summing up all momentum depositions we end up with transversally isotropic flow. The situation 
is less clear, however, in non-central collisions. The asymmetry of sizes between the direction in the 
reaction plane and that out of the reaction plane may be reflected in an asymmetry of the collective 
effect resulting from momentum deposition from hard partons. Such an effect is examined here. 

3. The effect of many jets 

The main idea here will be that the jets induce some kind of streams when they deposit mo- 
mentum into the bulk. In the literature, they are known as diffusion wakes. It is important that even 
if the leading partons are fully stopped, the streams seem to continue and carry momentum flfc)]. 
Originally, they flow in all transverse directions isotropically. Let us imagine, however, that two 
such streams come together from exactly opposite directions. Their momenta would cancel and 
energy would be deposited into the place of the merger. Thus the flow-generating effect would be 
smaller than just from a simple addition of two streams. In a very handwaving and cartoon-like 
way one could argue that in non-central collisions there is better chance of this to happen if the 
jets fly in the direction perpendicular to the reaction plane. The situation is sketched in Figure [[]. 
The fireball is elongated out of the reaction plane. Thus two jet-generated streams of bulk matter 
may have enough space to pass each other and not cancel, if they fly in the reaction plane. Then 
they both lead to observable asymmetry in the azimuthal hadron distribution. Streams which have 
their direction perpendicular to the reaction plane have less space available for passing by and the 
probability that they will merge and their momenta will cancel is higher. 

This reasoning would suggest, that the effect of having many hard partons inducing flow in 
the fireball and the individual flows merging would lead to a positive net contribution to the elliptic 
flow parameter vi- The proper way to test this conjecture would be a hydrodynamic simulation 
with the jets feeding the flow, i.e. technically source terms for in the hydrodynamic equations 
would be introduced. At the moment, such a simulation is technically too complicated with very 
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large number of jets. I therefore choose a simpler way and construct a toy model to represent the 
situation under study. 

4. The toy model 

The streams within the fluid are represented by blobs of matter. They all fly with velocities 
of 0.999c in various directions. Below it is described how the directions and initial positions are 
chosen. When two blobs meet, they merge into one which has larger mass and the momentum such 
that energy and momentum are conserved. In this way the merger of two streams (diffusion wakes) 
is represented. In the end, when there are no more mergers, the blobs evaporate pions according to 
a thermal distribution with a temperature of 170 MeV until all their energy is used up. 

Here we make a calculation for the LHC. 

The blobs carry momentum according to the distribution of hard partons in transverse energy 
and pseudorapidity, which is parametrised as 

^^•■O^p'Wv (4,, 
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The normalisation here is for single nucleon-nucleon collision. The differential cross section in 
Ej is integrated over the pseudorapidity interval [-2.5,2.5]. It is assumed that the distributions of 
jets in Ej and rj factorize in this interval. These parametrisations have been obtained from fits to 



MC results published in Ql lp - There, they were shown for Ej above 20 GeV. Here they will be 
extrapolated to lower Et's. As a result, the multiplicity of jets will be slightly overestimated. A 
more realistic parametrisation is being worked out. 

The azimuthal angle is distributed isotropically . The mass of the blobs is determined from the 
fixed velocity of 0.999c and the generated momentum. 

The total average number of blobs in one collision is given by the number of (mini)jets pro- 



duced. This can be calculated with the help of parametrisations (4.1) and (4.2). We first define the 



cross-section for a production of jets with transverse energy bigger than E, 



m 



a(E m )= T^dEr, (4.3) 
Je„, dE T 

and then determine the number of jets produced in the non-central nuclear collision of nuclei with 
mass number A as 

W) - pk^Ml (4 .4, 
In the last equation we have introduced the impact parameter b and the nuclear overlap function 

T M (b)= [ T A (r)T A (f-b)d 2 r, (4.5) 

./overlap 

which is defined with the help of the nuclear thickness function 



T A (r) = 2p jR 2 A - r 2. (4.6) 
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centrality 


b [fm] 




0-10% 


2.8 


14174 


10-30% 


5.9 


7820 


30-50% 


8.4 


3330 


50-75% 


10.5 


933 


75-100% 


12.4 


56 



Table 1: The centrality classes which were used in the simulation. For each class the impact parameter b 
and the average multiplicity of minijets with Ej > 4 GeV in the simulated interval Nj is presented. 

Integration in eq. ( |Q| ) goes over the whole overlap region of the two nuclei. In eq. (|4~6|), p 
is the nuclear density; it is assumed here that the density is uniform. Note also the factor K = 
10 fnr/jiib which converts the microbarns from the cross section parametrisation into fermis. 
The original transverse positions of the blobs are generated from the distribution 

p(r)=T A (r)T A (r-b), (4.7) 

which corresponds to scaling with the number of binary collisions. In the longitudinal direction 
and time, the space-time point of their appearance is determined as 

.x- = r u Q (4.8) 
x 3 = r m , (4.9) 

where is the four- velocity of the blob and To 1S the formation time of the stream, which is chosen 
here 0.6 fm/c in accord with the fast thermalisation conjecture [JIJ]. 

The size of the blobs can be varied in order to simulate various transverse sizes of the streams. 
It is decisive for the collisions because it actually gives the size of the cross section for blob-blob 
mergers. 

5. Results 

In the simulation from which results are presented the choice of parameters was E m = 4 GeV 
and the radius of the blobs 2.5 fm. Results were simulated in five centrality classes, which are 
summarised in Table |T} As mentioned and argued in the previous section, the jet multiplicities are 
most probably overestimated and better parametrisation of their production cross section is being 
elaborated. For the time being, let us proceed with the available parametrisation. 

In Figure ||]the azimuthal distribution of pions produced in this model is shown. In total, 10,000 
events were simulated. Histograms are shown for centrality classes 10-30% and 30-50%, where 
the azimuthal anisotropy of produced pion distribution is best seen. In both more and less central 
collisions the histogram is flatter or the statistics is worse. Though numerically the effect is small 
and gives V2 of just slightly less than 1% in this simulation, it is obvious from the histograms that 
the production asymmetry is correlated with the reaction plane. Note that the angle corresponds 
to the direction of the impact parameter. In other words, the angular distribution of pions produced 
in this model is given by the initial geometry and cannot be tagged as non-flow vi contribution. 
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Figure 2: The azimuthal distributions of produced pions for the centrality classes 10-30% (left) and 30-50% 
(right). 

6. Conclusions 

The numerical results presented here were obtained with a very simple toy model. Never- 
theless, the message that they carry is that there could be a contribution to the elliptic flow from 
feeding in momentum from the momentum loss of many hard partons which are expected at the 
LHC. It may be important to understand this contribution properly — together with all other effects 
determining the azimuthal asymmetry of the flow — if quantitative conclusions about the transport 
coefficients are to be drawn from the measurements. 

While we currently work on an improvement of the initial parametrisation of hard parton 
production, it goes far beyond the capability of the presented toy model to simulate proper response 
of the medium to momentum deposition from many hard partons. Interesting progress is being 
made in exploration of this problem, however. As argued here, it would be interesting to see a 
result of a proper hydrodynamic simulation of this situation. 
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